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FOREWORD 


This  Research  Note  reports  the  findings  of  The  llydroloqic  Fngineerinn 
Center  on  appropriate  methodoloqies  for  calculatinq  and  routinq  floods  re¬ 
sulting  from  suddenly-breached  dams. 

This  study  was  prepared  for  the  (J.S.  Army  Engineer  Waterways  Experiment 
Station,  Vickshuro,  Miss,  with  fundino  provided  bv  the  Defense  Nuclear 
Agency  under  subtask  L10HAXSX337,  "Above  Ground  Structures,"  work  unit  07, 


"Damage  of  Dams,"  and  by  the  Office,  Chief  of  Engineers  under  DA  project 
4A762710AT40,  task  A1 ,  work  unit  006. 

The  material  contained  herein  is  offered  for  information  purposes  onlv 
and  should  not  be  construed  as  Corps  of  Enoineers  policy  or  as  beinn  re¬ 
commended  guidance  for  field  offices  of  the  Corps  of  Engineers. 
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RFSFARCH  NOTE  NO.  5 

GUIDELINES  FOR  CALCULATING  AND  ROUTING  A  DAM-BREAK  FLOOn 


1.  Introduction .  Planning  and  design  requirements  for  a  wide  ranqe  of  pro¬ 
jects,  such  as  emergency  preparedness  and  siting  of  nuclear  power  plants, 
have  generated  widespread  interest  in  dam  break  floods.  Much  academic  re¬ 
search  and  some  laboratory  research  have  been  accomplished  on  this  topic. 
Generalized  analytic  techniques  for  calculating  and  routina  such  floods, 
particularly  in  non-prismatic  valleys,  have  not  been  readily  available. 
Furthermore,  prototype  verification  data  are  almost  non-existent.  This  re¬ 
port  describes  procedures  necessary  to  calculate  and  route  a  dam  break  flood 
usinn  an  existing  generalized  unsteady  open  channel  flow  model.  The  recent 
Teton  Dam  event  was  reconstituted  to  test  the  model's  performance  on  such  a 
highly  dynamic  wave.  The  procedures  outlined  herein  relate,  primarily,  to 
partial  breaches.  Some  deficiencies  in  the  model  were  identified  which  will 
require  some  further  research  and  programming  to  improve  the  applicability 
of  the  program  to  dam  break  flood  events. 

2.  Summary .  The  special  projects  nemo  cited  as  reference  (a)  established 
four  objectives  for  this  study.  The  first  two,  a)  level  of  accuracy  of 
existing  techniaucs  and  b)  sensitivity  of  calculated  results  to  n-values  and 
breach  size,  are  summarized  below  and  presented  in  detail  in  Appendix  A.  The 
third  objective,  c)  description  of  physical  phenomena  controllino  depth  and 
travel  time  and  a  discussion  of  pertinent  field  data,  is  presented  in  the 
body  of  this  report.  The  fourth  objective,  d)  documentation  of  the  method¬ 
ology,  is  included  in  Appendix  B.  Computer  programs  utilized  in  the  meth¬ 
odology,  references  (b)  and  (c),  na.v  he  obtained  from  The  Hydrologic 
Engineering  Center. 

The  computer  program  of  reference  (c)  was  applied  to  the  Teton  Dam  data  set 
to  demonstrate  the  level  of  accuracy  one  mioht  expect  in  such  analyses. 

The  results  are  shown  on  pages  A-?fi  through  A-?8  of  Appendix  A  and,  in 
general,  appear  reasonable.  This  test  case  demonstrates  the  usefulness  of 
a  generalized  computer  program  because  the  methods  proposed  in  references 
(d)  and  (e)  were  not  applicable  to  the  Teton  data  set  for  reasons  given  in 
paragraph  10. 


Regarding  sensitivity  to  breach  size,  pages  A-22  and  A-23  show  the  two 
breach  sizes  considered.  The  breach  that  developed  at  Teton  was  estimated, 
by  others,  to  be  40  percent  of  the  dam  embankment.  Geometric  data  were  not 
available  to  verify  this,  therefore,  our  best  estimate  of  the  final  Teton 
breach  geometry,  page  A-22,  is  based  on  photographs.  The  breach  shown  on 
page  A-23  has  the  same  side  slope  as  that  on  page  A-22,  0.6  on  1  but  It 
has  zero  bottom  width.  This  seemed  a  likely  intermediateconditfon,  but 
no  field  data  were  available  at  the  time  of  this  study  to  establish  an 
observed  intermediate  condition. 

The  calculated  outflow  is  shown  on  page  A-24.  The  hydrograph  labeled 
"trapezoidal  breach"  assumed  the  40  percent  breach  size,  page  A-22, 
developed  instantaneously.  The  hydrograph  labeled  "triangular  breach" 
was  determined  in  a  similar  manner  for  the  TO*?  breach  size.  The  third 
hydrograph  on  page  A-?4  was  calculated  for  the  trapezoidal  breach  (labeled 
40?,  breach  size  on  page  A-'72),  but  an  observed  reservoir  drawdown  curve 
at  the  dam,  page  A-20,  was  used  which  implies  a  gradual  development  of 
the  breach  rather  than  instantaneous  failure.  The  last  approach  was  con¬ 
sidered  best  in  estimating  the  discharge  hydrographfrom  Teton  reservoir 
given  the  data  set  and  analytical  techniques  available  to  us. 

The  sensitivity  of  calculated  outflows  to  breach  size  and  rate  of  development 
is  illustrated  on  page  A-24.  It  is  summarized  in  the  following  table  together 
with  pertinent  elevation  data  for  an  n  value  of  0.04. 


Table  1:  Sensitivity  to  Breach  Size  and  Rate  of  Development 


Final  Breach  Size 
%  of  Total  Dam 

Rate  of 
Development 

[calculated  Peak 
Hater  Discharge 
at  Dam  Axis 

106  CFS  (1) 

Calculated  Peak  Elevations 

MSU 

At  Dam  Axis 

& 1 

171 

T’r) 

ITH  (3) 

ittuyH 

TO* 

instanta- 

2.4 

5151 

5015 

4933 

neous 

40? 

instanta- 

3.4 

5175 

5020 

4935 

neous 

(1)  Multiply  by  0.0?032  to  get  Cubic  deters  Per  Second 

(2)  Actual  rate  of  development  was  unknown  so  the  observed  reservoir 
drawdown  curve,  page  A-20,  was  used  to  approximate  outflow 
conditions. 


(3)  The  actual  peak  discharge,  as  estimated  bv  personnel  of  the  Walla 
Jalla  District,  Armv  Corps  of  Engineers  from  observed  data  in  the 
Teton  Canyon  three  miles  downstream  from  the  dam,  was  2,300,000  rfs. 


From  these  results  it  is  apparent  that  neither  the  size  of  breachs  tested 
nor  the  rates  of  failure  assumed  were  very  significant  in  predicting  peak 
elevations  five  miles  downstream  from  the  dam. 

The  calculated  peak  flood  elevations,  near  the  dam,  were  very  sensitive 
to  n-values.  Increasing  n  from  .03  to  .06  raised  the  peak  flood  elevation 
25  feet  at  the  dam,  as  Illustrated  on  page  A-8.  TabV  At  5  miles 
downstream  the  calculated  difference  was  only  8  feet.  Differences  continued 
to  diminish  with  distance. 

Calculated  Travel  Times  are  shown  on  page  A-29.  They  correspond  to  the 
discharge  hydroqraph  labeled  "simulated  from  observed  data"  on  paqe  A-24 
and  n-values  of  0.04. 

Searching  for  a  simolified  approach  in  place  of  references  (d)  and  (e)  led  to 
a  trial  application  of  the  Modified  Puls  routing  technique.  The  hydroqraph 
labeled  "simulated  from  observed  data"  on  page  A-24  was  routed  and  a  water 
surface  profile  calculated  for  the  resulting  peak  discharges.  A  comparison 
of  the  results  with  the  observed  elevations  and  the  peak  elevations  computed 
with  the  full  equations  is  shown  on  pages  A- 30  through  A-32.  Additional 
investigation  is  needed  to  establish  the  range  of  applicability  of  tins  * 
method. 


3.  Physical  Phenomena  and  Field  Data.  Analysis  of  the  dam-break  flood 
involves  understanding  the  physical  processes  before  aoplyinq  analytical 
techniques  which  approximate  those  physical  processes.  Three  distinctly 
different  processes  are  involved:  the  process  of  structural  failure 
causing  the  breach  to  develop;  the  process  of  setting  water  into  motion 
in  a  reservoir;  and  the  process  of  flood  wave  attenuation. 

The  size,  shape  and  rate  of  breach  development  are  primarily  responsible 
for  the  peak  rate  of  outflow  from  the  reservoir.  Yet,  of  the  three  physical 
processes,  this  one  is  the  most  difficult  to  quantify.  With  the  exception 
of  man-made  breaches,  it  is  difficult  to  visualize  the  instantaneous  develop¬ 
ment  of  a  breach.  Some  have  occurred,  however.  The  St.  Frances  Dam,  a 
high  head  concrete  gravity  structure,  apparently  suffered  an  abutment 
failure  which  resulted  in  virtually  the  instantaneous  failure  of  the  entire 
structure.  The  Johnstown  flood  of  1889  was  caused  by  the  complete  failure 
of  an  earth  fill  dam.  Reports  Indicate  that  less  than  half  an  hour  was 
required  for  overtopping  flow  to  breach  the  structure.  The  recent  Teton 
failure,  a  full  depth-partial  width  breach  of  an  earth  fill  dam,  is  estimated 
to  have  developed  in  less  than  two  hours.  Since  natural  failure  of  a  major 
structure  Is  so  Improbable,  establishing  a  mode  of  failure  requires  a  policy 
decision  rather  than  an  analytical  technique.  In  qeneral ,  instantaneous 
failure  of  the  entire  structure  produces  the  largest  flood  wave. 
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The  second  physical  process  results  from  the  depth  of  water  above  the  breach 
invert.  That  Is,  a  reservoir  has  a  total  enerqy  head  equal  to  the  elevation 
of  the  water  surface.  If  the  dam  Is  breached,  the  force  of  qravity  will  set 
water  into  motion.  The  effect  will  propogate,  as  a  negative  wave,  to  the 
upstream  end  of  the  reservoir  at  a  velocity  equal  to  yfqy  where  g  Is 
acceleration  of  gravity  and  y  vs  water  depth.  Because  of  the  great  depth 
in  a  reservoir,  very  little  frictional  resistance  is  mobilized  during  the 
passage  of  this  negative  wave.  As  a  result,  water  gains  specific  energy 
rapidly  as  It  moves  toward  the  breach.  In  instantaneous  breach  development, 
the  peak  outflow  will  occur  within  a  minute  or  two  after  breaching. 

Whereas  the  total  energy  head  setting  the  water  into  motion  is  the  specific 
energy  (i.e.,  the  initial  water  depth)  above  the  breach  invert,  the  energy 
which  must  be  dissipated  in  the  downstream  channel  is  equal  to  the  specific 
energy  from  the  downstream  channel  invert  to  the  initial  pool  elevation. 

The  fact  that  the  water  surface  elevation  drops  down  rapidly  at  the  dam 
axis  does  not  reflect  a  corresponding  loss  in  energy  head.  When  flow 
begins,  that  specific  energy  above  the  breach  invert  is  transformed  into 
three  components:  a  pressure  head,  a  kinetic  energy  head  and  an  inertia 
head.  (The  relative  size  of  each  of  these  energy  head  components  is 
discussed  more  fully  in  sections  4  and  5.)  Friction  loss  is  relatively 
small  and  may  be  neglected  unless  the  reservoir  bottom  is  extremely  rough 
(more  than  5"  or  10"  of  the  water  depth). 

The  third  physical  process,  flood  wave  attenuation,  involves  energy 
dissipation  and  valley  storage.  As  the  flood  wave  moves  downstream,  the 
peak  discharge  tends  to  decrease,  the  base  of  the  flood  wave  will  become 
longer  and  the  wave  velocity  will  decrease.  Hear  the  dam,  enerqy  dissipation 
is  primarily  responsible  for  behavior  of  the  flood  wave.  However,  valley 
storage  soon  becomes  the  primary  factor  in  flood  wave  attenuation.  The 
key  to  the  transition  from  energy  dissipation  to  valley  storage  control  is 
the  rate  at  which  the  slope  of  the  total  energy  gradient,  a  line  which  must 
intersect  the  initial  pool  elevation  at  the  dam,  is  reduced  to  that  of  a 
major  rainfall  flood  in  the  downstream  valley.  It  seems  obvious  that  the 
total  energy  at  any  cross  section  in  the  valley  should  not  exceed  the  initial 
reservoir  elevation,  and  yet  some  analytical  techniques  occasionally  violate 
that  constraint.  It  is  good  policy  to  always  check  the  total  energy,  as 
well  as  the  water  volume,  in  a  calculated  flood  wave. 

The  rate  of  energy  dissipation  is  qoverned  primarily  by  friction  loss. 

Minor  losses  from  bends  and  contractions-expansions  are  often  included  in 
the  n-values. 

The  volume  of  water  in  the  reservoir  is  the  final  piece  of  field  data  required. 
This  volume  strongly  Influences  the  peak  elevations  at  downstream  points. 


4.  Energy  Components  and  Peak  Outflow  from  Complete,  Instantaneous  Breaches. 
It  is  useful  to  develop  the  relative  size  for  each  energy  component  in  the 
flow  at  the  dam  axis  and  to  compare  all  of  them  to  the  more  conmon  case  of 
steady  state  critical  flow  at  a  contraction. 

By  assuming  a  rectangular  cross  section,  zero  bottom  slope  and  instantaneous 
removal  of  the  entire  dam,  Saint-Venant  developed  an  analytical  solution  for 
the  elevation  of  the  free  surface,  reference  (f).  page  755.  Utilizing  that 
equation,  the  depth  of  flow  at  the  dam  axis  was  determined,  by  Saint-Venant 

and  others,  to  be  ^  Y0  where  Y0  is  the  original  water  depth  at  the  dam.  _ 

Also,  the  velocity  corresponding  to  the  peak  outflow  was  shown  to  be  \  VqY0. 
Combining  these  relationships  leads  to  the  equation  for  peak  discharge 


q 


8 

max  ’  Tf 


(1) 


Yq  is  the  initial  water  depth  at  the  dam 
q  is  acceleration  of  gravity 
g  is  peak  water  discharge  in  cfs/ft 

m3  X 


Since  this  equation  was  developed  for  a  rectanqular  section,  the  total 

discharge  may  be  calculated  by  multiplying  q  by  the  width. 

m3  x 

Using  the  relationships  referenced  above,  the  velocity  head  (i.e.,  the 
kinetic  energy  head  component  of  the  specific  energy  head)  was  calculated 
to  be  y-  Y„.  Since,  in  the  absence  of  friction  and  other  losses,  inertia 
is  the  only  remaining  term  in  the  basic,  unsteady  flow  equations  of  Saint- 
Venant,  it  may  be  calculated  as  follows. 
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These  components  are  shown  in  Figure  1  along  with  the  energy  components 
for  critical,  steady  state  flow. 

This  figure  shows  that  in  the  dam  break  flood  analysis,  as  well  as  steady 
state  critical  flow  at  a  contraction,  the  velocity  head  is  half  the  pressure 
head.  However,  the  inertia  head  comnonent  is  zero  in  Figure  la  because  flow 
is  steady  state. 


r. 


(a)Critical  flow  at  a  contraction  (h)Critical  flow  from  a  breached  dam 

(steady  flow  conditions)  Unsteady  *lo,-<  c'w<it1nnsl 

Figure  1.  Components  of  Specific  Energy  Head. 


The  drawdown  in  water  surface  elevation  to  4"  Y0  at  the  dam  axis,  Fiqure  lb, 
does  not  reflect  a  corresponding  energy  loss.  Experimental  results  obtained 
by  Schoklitsch,  reproduced  on  page  755  of  reference  If,  show  relatively 
little  friction  loss  in  flow  approachinq  the  dam  axis.  As  might  be  expected, 
the  model  results  showed  friction  to  be  very  significant  downstream.  Tests 
reported  hy  WES  in  reference  (g)showed  no  Impact  from  friction  loss  at  the 
dam  axis.  However,  the  WES  flume  sloped  at  0.005  ft/ft,  whereas  the  flume 
in  Schokl itsch1 s  experiment  had  zero  bottom  slope. 

The  significance  of  this  point  is  that  all  three  enerqy  components,  pressure 
head,  kinetic  enerqy  head  and  inertia  head,  are  significant  in  complete, 
instantaneous  breachinqs.  Consequently,  investiqators  encouraqe  the  use  of 
the  complete  routing  equations,  often  referred  to  as  the  Saint- Venant 
equations.  Simplifications  of  the  complete  equations,  such  as  Musklngham, 
Tatum,  Straddle-Staqqer  and  'V>di f i ed  Puls,  are  not  recommended  because  the 
empirical  coefficients  would  invariably  be  developed  from  rainfall  floods 
and  would  reflect  different  values  of  energy  components  relative  to  Y0. 

5.  Instantaneous.  Partial  Breaches.  Partial  breaches  are  classified, 
according  to  hydraulic  performance,  as  full  depth-partial  width,  partial 
depth-full  width  or  partial  depth-partial  width.  A  separate  equation  has 
been  developed  for  calculating  the  peak  water  discharge  for  each  class, 
page  25  of  reference  (g). 
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An  empirical  equation  for  partial  depth-partial  width  breaches  was  reported 
in  references  (g)  and  (h). 
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For  breach  sizes  in  the  following  range. 
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Since  the  discharge  equations  for  partial  breaches  are  similar,  in  form, 
to  that  for  a  full  breach  (1),  the  total  specific  energy  has  the  same  three 
basic  components.  However,  their  size,  relative  to  initial  water  depth, 
is  considerably  different  from  that  shown  in  Figure  1.  There  is  no  analytical 
solution  for  partial  breaches,  therefore,  experimental  results,  presented  in 
reference(g) ,  were  used  to  calculate  the  individual  energy  head  components. 

The  following  table  presents  experimental  results  for  full  depth  breaches 
ranqing  in  width  from  10 %  to  100%  of  the  flume  width  in  columns  1,  2  and  3. 
Fractions  of  initial  water  depth,  calculated  with  equation  3,  are  shown  in 
columns  4  and  5.  A  sample  of  the  calculations  is  presented  in  the  paragraph 
following  the  table.  This  sample  calculation  utilizes  equation  3  and  a  100% 
breach  size  (i.e.,  full  breach)  to  demonstrate  that  the  relative  value  of  each 
energy  component  is  the  same  as  the  respective  value  produced  by  equation  1, 
the  analytical,  full  breach  equation,  when  equation  3  Is  carried  to  Its 
upper  limit. 


Table  ?:  Relative  Size  of  Enerqy  Components  in  Partial  Width  Breaches^ 


Test 

No 

Breach 

Size 

(V 

Pressure 

Head 
<  Of  Yq 

Velocity 

Head  (2) 

S  of  Y0 

Inertia 
Head  (3) 
%  of  Y 

0 

(1) 

(2) 

(3) 

(4) 

(5) 

1.1 

rull 

44 

?2 

34 

2.1 

60 

70 

12 

18 

3.1 

30 

8? 

12 

6 

4.1 

15 

89 

(4), 

-  - 

5.1 

10 

94 

(4) 

-  - 

Notes:  1.  Values  in  columns  1  and  2  are  from  Table  A,  paqe  8,  reference 
(a)and  values  in  column  3  are  from  experimental  results  from 
Tables  1  through  5,  Station  200,  reference  (p). 

2.  Velocity  head  is  calculated  with  equations  8  and  9,  following. 

3.  Inertia  head  is  Y0  -  (pressure  head  +  velocity  head). 

4.  Calculated  values  exceeded  100  percent  of  Y0,  which  probably 
reflects  scatter  in  experimental  results. 
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0  from  equation  (3) 
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For  the  full  breach,  b  =  1.0B  and  y  =  ^  y 
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This  agrees  with  section  4  and  shows  the  procedure  followed  In  completing 
Table  2.  The  inertia  head,  column  5  in  Table  2,  was  calculated  assuming 
zero  energy  loss  upstream  from  the  dam. 


h.- 
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Recause  of  the  decrease  in  relative  significance  of  inertia  head  and  even 
velocity  head,  it  is  satisfactory  to  apply  simplifications  of  the  full 
Saint-Venant  equations  to  oartially  breached  dams. 

6.  Attenuation  of  the  Flood  Have.  As  a  flood  wave  moves  downstream, 
friction  and  other  losses  chanqe  the  relative  size  of  the  three  energy 
components.  Even  floods  from  fully  breached  dams  eventually  take  on  the 
characteristics  of  i  rainfall  flood  and  may  be  routed  with  a  simplified 
routing  method  such  as  Modified  Puls.  Major  areas  of  uncertainty  are  1) 
how  much  distance  is  required  for  this  transition,  2)  how  does  this  distance 
vary  when  considering  partial  breaches  and  3)  what  is  the  maximum  breach 
size  to  consider  as  a  partial  breach. 

7.  Proposed  Analytical  Technique.  The  guidelines  presented  in  Appendix 
R  of  this  report  are  developed  for  the  computer  program  "Gradual ly  Varied 
Unsteady  Flow  Profiles".  It  is  a  solution  of  the  basic  Saint-Venant 
equations  for  unsteady  flow  and  may  be  used  to  calculate  the  outflow 
hydrograph  through  any  size  or  shape  of  breach,  as  well  as  to  route  that 
hydroqraph  downstream  and  provide  water  discharge  and  water  surface  elevation 


9 


hydroqraphs  at  any  number  of  computation  points  up  to  45.  The  maximum 
discharge,  maximum  elevation  and  maximum  flow  velocity  are  summarized  for 
each  computation  point. 

Sufficient  information  is  printed  out  so  the  time  of  arrival,  time  of  peak 
and  duration  of  the  flood  may  be  plotted. 

This  computer  program  accounts  for  the  movement  of  the  negative  wave  through 
the  reservoir,  for  the  tail water  submergence  at  the  dam,  for  the  three 
components  of  energy  presented  earlier,  for  friction  loss  and  for  storaqe 
in  the  reservoir  and  the  downstream  valley. 

Cross  sections  need  not  be  rectangular  or  Drismatic.  A  companion  program, 
"Geometric  Elements  from  Cross  Section  Coordinates",  is  available  to 
transjVm  complex  cross  sections  into  the  required  qeometric  data  set  for 
the  routing  proqram. 

These  computer  programs  are  generalized.  That  is,  they  are  sufficiently 
flexible  and  adaptable  to  be  used  without  code  changes.  They  are  portable 
from  one  computer  to  another  and  documentation  is  available,  from  The 
Hvdroloaic  Engineering  Center. 

8.  Program  Limitations. 

a.  Routing  with  the  Gradually  Varied  Unsteady  Flow  Profiles  computer 
program  requires  a  large  high  speed  computer  (50,000,  fiO-blt  words)  and 
personnel  who  are  experienced  in  applying  mathematical  models. 

b.  Any  breach  size  nay  be  modeled,  but  the  proqram  assumes 
instantaneous  development. 

c.  All  channels  must  be  wet  initially.  That  is,  computations  cannot 
be  made  if  any  portion  of  the  model  is  dry.  This  is  overcome  by  prescribing 
a  base  flow;  however,  the  computer  program  has  difficulty  in  establishing 
this  profile. 

d.  ’’ovement  of  the  negative  wave  through  the  reservoir  causes  no 
computational  problem  until  it  reaches  the  upstream  end  of  the  reservoir. 
Computation  nodes  tend  to  go  dry  and  abort  the  computer  run. 

e.  The  analysis  of  multiple  failures  would  require  manual  intervention 
to  stop  and  restart  the  calculation  process  as  each  new  structure  is  brought 
into  the  system. 

f.  The  program  assumes  a  horizontal  water  surface  transverse  to  the 
flow,  whereas  a  great  deal  of  transverse  slope  can  exist  in  the  actual 
prototype  situation. 


y.  Proposeu  Areas  of  Research.  All  of  the  program  limitations  were 
circumvented  in  analyzing  the  Teton  Data  Set.  The  trade-off,  however,  was 
analysis  time.  Seven  weeks  were  required  to  set  up  the  data,  debug  it  and 
perform  the  analysis.  The  two  tasks  requiring  the  most  time,  probably  75?', 
were  establishing  initial  base  flow  conditions  for  the  model  (8c)  and 
stabilizing  the  computations  when  the  negative  surge  reached  the  upstream 
boundary  (Bd) .  Both  of  these  problem  areas  can  be  overcome  by  additional 
programming.  The  improvements  would  reduce  analysis  time  to  four  or 
five  weeks. 


Instantaneous  breach  development,  8b,  could  be  replaced  by  equations  which 
let  progressive  development  take  place.  In  the  absence  of  a  theory,  the  rate 
of  development  would  have  to  be  prescribed  with  input  data. 


Developing  the  capability  to  handle  multiple  dam  failures  (3e),  especially 
in  tandem,  will  he  a  major  modification. 


This  analytical  technique  is  a  one-dimensional  model  and  will  always  have  a 
Horizontal  water  surface  transverse  to  the  flow.  At  present,  two-dimensional 
modeling  is  not  feasible. 

ID.  A1 ternate  Analytical  Procedures .  Alternate  analytical  procedures  were 
proposed  in  references  (d)  and  Je) ■  Neither  were  applicable  to  the  Teton 
Data  Set. 

The  dimensionless  curves  were  developed  from  numerical  solution  of  the  St. 
Venant  equations  and  include  special  treatment  of  the  wave  front  as  it  moves 
along  a  dry  channel.  By  knowing  reservoir  volume,  valley  cross  section  at 
the  dam,  initial  reservoir  elevation,  stream  slope  and  stream  roughness, 
the  curves  will  pro  vice  three  properties  of  the  flood  wave: 

1.  Time  of  arrival  at  downstream  points 

2.  Maximum  depth  profile  in  the  downstream  channel 

3.  Time  of  maximum  depth  at  downstream  points. 

The  curves  extend  for  distances  ranging  up  to  fifteen  times  the  reservoir 
length.  The  outflow  hydrograph  at  the  dam  is  not  needed  to  use  these 
curves.  It  was  assumed,  in  developing  the  curves,  that  the  entire  dam 

is  breached  instantaneously  and  that  the  valley  is  prismatic.  Neither 
conoition  was  satisfied  by  the  Teton  case. 

Tne  procedure  in  reference  (e)was  developed  for  smaller  structures  and  the 
Teton  Data  Set  was  completely  beyond  the  range  of  nomographs  and  curves 
presented  there.  In  any  case,  the  procedure  does  not  route  the  flood  wave 
downstream.  Only  the  outflow  discharge  hydrograph  is  calculated  at  the 
uan  axis.  The  procedure  can  handle  a  wide  range  of  breach  sizes,  but  it 
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is  designed  with  partial  breaches  in  mind.  It  has  the  advantage  of  tail 
water  correction,  which  is  essential  when  breaching  of  a  low  dam  coincides 
with  a  high  flow  condition  In  the  stream.  The  procedure  is  well  documented 
and  Is  easily  applied. 

A  possible  alternative  approach  for  partial  breaches  is  the  Modified  Puls 
routing  technique.  Preliminary  work  with  this  technique  produced  the  re¬ 
sults  shown  on  pages  A- 30  through  A-32  for  the  Teton  Data  Set.  A  Manning 
n  value  of  0.04  was  used;  further  details  are  given  In  Paragraph  5,  Appendix 
A.  The  advantage  of  this  technique  is  that  readily  available  and  easily 
applied  computer  programs  (e.g.,  HEC-1  and  HEC-2)  can  be  utilized;  total 
analysis  time  would  probably  be  reduced  to  two  to  three  weeks. 


The  disadvantage  is  that  the  range  of  application  is  limited  whereas  the 
technique  presented  in  Paragraph  7  is  generally  applicable. 

Additional  research  is  needed  to  define  the  range  of  applicability  of  the 
Modified  Puls  technique.  The  present  hypothesis  is  that  the  size  of  the 
inertia  component.  Table  2,  would  provide  a  suitable  parameter  for  defining 
that  range. 

This  research  would  not  require  additional  physical  modeling.  Studies  re¬ 
ported  in  references  (g)  and  (h)  offer  test  data  for  numerical  studies. 

Other  numerical  experiments  could  be  performed  by  using  results  from  anal¬ 
yzing  variation s  of  the  Teton  Data  Set  with  the  complete  equations.  These 
results  could  be  obtained  while  pursuing  any  of  the  areas  of  research  pro¬ 
posed  in  Section  9. 

Computer  programs  which  utilize  the  Modified  Puls  routing  technique  are 
available  and  are  presently  developed  to  a  higher  degree  of  serviceability 
than  programs  solving  the  full  equations.  Water  surface  profile  computations 
will  be  required  In  conjunction  with  the  Modified  Puls  routing  to  produce  a 
water  surface  profile.  These  computations  are  computerized  also.  No  major 
computer  program  development  would  be  required.  The  appropriate  existing 
computer  programs,  HEC-1  and  HEC-2,  are  well  documented. 
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appendix  a 


UNSTEADY  FLOW  ANALYSIS 
TETON  DAM  FAILURE 


1.  INTRODUCTION 

This  phase  of  the  study  which  calculates  and  routes  the  flood  resulting 
from  a  suddenly  breached  dam,  consists  of  a  computer  solution  in  conjunction 
with  the  Teton  Dam  failure  as  defined  in  objectives  a  and  b  of  Special 
Projects  Memo  No.  473.  The  analysis  utilizes  the  unsteady  flow  computer 
program  to  determine  water  surface  elevations  resulting  from  various  breach 
sizes  and  n  values.  The  level  of  accuracy  was  determined  by  comparing  avail¬ 
able  flood  data  (particularly  high  water  marks)  from  the  June  5,  1976  dam 
failure  with  calculated  results. 

The  analysis  used  data  generally  available  to  field  personnel 
engaged  in  the  study  of  the  impact  of  a  dam  break  flood  such  as  topo¬ 
graphic  maps,  aerial  photography,  dam  description,  gaged  stream  flows 
and  reservoir  elevation-capacity  curves.  The  primary  area  of  study 
included  about  30  miles  of  flood  plain  downstream  of  the  dam. 

2.  GENERAL  DESCRIPTION 

The  Teton  Dam  is  located  on  the  Teton  River  in  southeastern  Idaho 
approximately  13  miles  north  and  east  of  the  city  of  Rexburg  (see  Plate  1). 
The  dam  was  designed  as  a  zoned  earthfill  embankment  with  a  crest  elevation 
of  5,332  feet  (mean  sea  level  datum)  and  a  maximum  height  of  305  feet 
(above  riverbed).  It  would  create  a  reservoir  of  288,250  acre-feet  when 
filled  to  an  elevation  of  5,320  feet.  (Plan  and  sections  of  the  dam  are 
shown  in  Plates  2  and  3.)  The  dam  is  located  in  a  narrow  steepwalled 
canyon.  The  channel  geometry  is  essentially  the  same  for  some  4  miles 
downstream  of  the  dam  whereupon  the  Teton  River  enters  upon  a  wide 
relatively  flat  flood  plain. 

The  total  reservoir  storage  just  prior  to  failure  on  June  5  was  about 
251,300  acre-feet  at  an  elevation  of  5,301.5  feet.  Measured  inflow  was 
3,580  cubic  feet  per  second  and  measured  outflow  before  any  leaks  developed 
was  940  cubic  feet  per  second. 

Although  water  stored  upstream  of  the  dam  caused  the  flooding,  other 
sources  contributed  significantly  to  the  floodflow  in  the  downstream  reaches 
However,  the  complexity  of  major  irrigation  diversions  and  numerous  return 
flows  precludes  an  accurate  accounting  of  all  sources.  According  to  the 
record  obtained  from  the  damaged  recorder  at  the  gaging  station  on  Henrys 
Fork  near  Rexburg  ( 


were  flowing  in  Henrys  Fork  at  the  time  the  leading  edge  of  the  flood 
wave  arrived  about  4:00  p.m.  on  June  5.  The  contribution  of  the  Snake 
River,  just  upstream  of  its  confluence  with  Henrys  Fork  was  estimated 
at  5,600  cubic  feet  per  second  around  1:00  p.m.  of  the  same  day. 

Table  1  presents  a  preliminary  tabulation  of  available  data  pertaining 
to  the  leading  edge  of  the  flood  wave  generated  by  the  dam  failure. 

Distances  shown  are  measured  from  the  dam  breach  along  the  estimated 
path  of  the  leading  edge  of  the  wave. 

Instantaneous  peak  discharges  were  determined  at  two  sites  along  the 
study  reach  (Table  2).  Indirect  measurements  based  on  field  surveys  and 
empirical  formulas  were  used  to  compute  the  peak  discharges  at  the  sites. 

The  indirect  methods  for  computing  peak  discharge  were  based  on  hydraulic 
equations  relating  discharge  to  the  water-surface  profile  and  the  geometry 
of  the  channel . 

3.  GEOMETRIC  MODELS 

The  channel  configurat ion  upstream  and  downstream  of  the  dam  was 
modeled  independently  according  to  procedures  contained  in  the  users 
manual,  generalized  computer  program  (723-G2-L745B) ,  Geometric  Elements 
from  Cross  Section  Coordinates  (GEDA),  Hydrologic  Engineering  tenter,  (J.S. 

Army  Corps' of  Engineers,  June  1976. 

The  upstream  geometric  model  (reservoir)  was  developed  essentially  from 
two  cross  sections  located  at  the  damsite  and  at  the  estimated  upstream 
limit  of  the  reservoir  (when  filled  to  capacity).  The  cross  section  data  at 
the  damsite  were  obtained  for  preproject  conditions  from  topographic  informa¬ 
tion  shown  on  Plates  2  and  3.  Information  for  the  other  section  was  obtained 
from  U.S.  Geological  Survey  topographic  maps  (scale  1:24,000).  Data  were 
coded  in  the  standard  HEC-21  format  and  input  into  the  GEDA  program.  Computed 
results  were  compared  with  respective  values  taken  from  the  reservoir  capacity 
curve  of  Plate  2  and,  where  necessary,  adjustments  made  to  the  upstream  cross 
section  until  computed  volumes  (at  specified  elevations)  plotted  reasonably 
close  to  the  published  curve  (see  Plate  4). 

The  downstream  model  (channel)  was  developed  in  a  similar  manner  although 
a  more  detailed  definition  of  the  channel  configuration  was  necessary.  Since 
the  actual  flooded  area  resultinq  from  the  dam  failure  included  several 
major  tributaries  (the  Snake  River,  Henrys  Fork  and  the  Teton  River),  there 
was  no  readily  available  information  on  storage  (volume)  within  the  study 
reach  with  which  to  calibrate  the  geometric  model.  Cross  sections  were 
taken  on  the  average  every  mile,  and  conveyance  limits  were  designated 
so  as  to  effectively  model  conveyance  for  the  large  flows  expected  in  the 
dam  break  analysis. 


*HEC-2  Water  Surface  Profiles,  Generalized  Computer  Program  (723-02A), 
Hydrologic  Engineering  Center,  U.S.  Army  Corps  of  Engineers,  October  1973. 
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TABLE  1 


ESTIMATED  TIMES 

FOR  LEADING  EDGE  OF  FLOOD  WAVE 
(BASED  ON  FIELD  OBSERVATIONS) 


Location 

Miles3 
downstream 
from  dam 

Date 

(month/day ) 

Arrival 

time 

(hours) 

Elapsed  time 
between  sites 
(minutes) 

Approximate 
mean  velocity 
between  sites 
(ft/min)(mi/hr) 

At  damsite 

6-5 

1157 

— 

— 

— 

In  Teton  Canyon 

3.0 

6-5 

8 

1  ,980 

23 

At  Teton 

8.8 

6-5 

25 

1  ,220 

14 

At  Sugar  City 

12.3 

6-5 

30 

620 

7 

At  Rexburg^ 

15.3 

6-5 

1340 

40 

5 

Henrys  Fork 
near  U.S.G.S. 
gaging  station, 
Henrys  Fork  near 
Rexburg 

22.6 

6-5 

1530 

110 

350 

4 

Menan  Bridge  on 
the  Snake  River 
immediately 
upstream  of  the 
Union  Pacific 

30.6 

6-5 

1800 

150 

280 

3 

distances  measured  along  the  estimated  path  of  the  leading  edge  of  the  flood  wave. 
bMaximum  depths  estimated  at  6-8  feet  within  the  city  of  Rexburg. 
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TABLE  2 


ESTIMATED  PEAK  DISCHARGES 
(BASED  ON  FIELD  OBSERVATIONS) 


Location 

Miles3 
downstream 
from  dam 

Date 

(month/day) 

Time 

(hours) 

Discharge 

(cubic  feet  per  second) 

In  Teton  Canyon 

3.0 

6-5 

b 

At  Teton 

8.8 

6-5 

c 

distances  measured  along  the  estimated  path  of  the  leading  edge  of  the  flood  wave, 
kpeak  probably  occurred  between  1230  and  1330  hours. 
cPeak  probably  occurred  between  1300  and  1400  hours. 
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The  area,  hydraulic  radius,  top  width  and  average  n  value  were 
calculated  at  each  cross  section  in  the  goemetric  models  for  the  eleva¬ 
tions  specified.  By  calculating  length-weighted  values,  the  preceding 
elements  were  modified  such  that  they  would  apply  to  uniformly  spaced 
computation  points  along  the  length  of  each  model  (reservoir  and  channel). 

4.  UNSTEADY  FLOW  MODELS 

The  unsteady  flow  data  models  of  the  reservoir  and  the  channel  were 
developed  according  to  procedures  contained  in  the  users  manual,  generalized 
computer  program  (723-G2-L2450) ,  Gradually  Varied  Unsteady  Flow  Profiles, 
Hydrologic  Engineering  Center,  U.S.  Army  Corps  of  Engineers,  January  1976. 

Reservoir  Model 


In  the  case  of  the  reservoir,  tables  of  geometric  elements  were  prepared 
at  21  uniformly  spaced  computation  points  (approximately  2,800  feet  apart) 
and  input  into  the  unsteady  flow  program.  Since  inflow  into  the  reservoir 
just  prior  to  failure  was  measured  (approximately  3,600  cubic  feet  per 
second),  the  upstream  boundary  condition  was  based  on  such.  A  normal 
depth  was  determined  for  the  measured  discharge  rate  at  the  upstream  cross 
section.  A  steady  flow  condition  was  assumed,  and  a  constant  elevation 
(determined  from  the  normal  depth  calculation)  was  specified  as  the  boundary 
condition.  Initial  values  (discharge  and  elevation)  were  specified  only 
at  the  most  upstream  and  most  downstream  computation  points  (nodes)  and 
a  linear  interpolation  scheme  within  the  program  utilized  to  determine 
values  at  the  remaining  nodes.  The  values  of  discharge  and  elevation  at 
the  upstream  node  are  the  same  as  those  utilized  in  the  computation  of  the 
elevation  hydrograph  at  the  upstream  boundary.  The  elevation  used  at  the 
downstream  node  (damsite)  was  5,301.5  feet  (the  observed  reservoir  water 
surface  elevation  at  time  of  failure)  and  the  discharge  was  estimated 
as  3,000  cubic  feet  per  second.  (Although  the  measured  discharge  was  940 
cubic  feet  per  second  before  any  leaks  occurred,  a  higher  discharge  rate 
was  used  because  of  significant  seepage  just  prior  to  failure).  In  the 
reconstruction  of  the  actual  discharge  hydrograph  at  the  damsite,  the 
observed  reservoir  water  surface  elevations  at  particular  points  in  time 
(after  failure)  were  input  initially  for  the  downstream  boundary  condi¬ 
tion  (see  Plate  5  for  a  plot  of  the  observed  data).  Since  the  observed 
data  were  not  adequate  with  respect  to  time  to  define  a  reasonable  eleva¬ 
tion  hydrograph,  a  greater  number  of  coordinate  points  were  coded  in  the 
final  analysis  (refer  to  Plate  5).  The  computed  discharge  hydrograph  at 
the  damsite,  based  on  the  preceding  conditions,  is  shown  in  Plate  6.  (Note 
that  the  computed  maximum  discharge  is  about  1,800,000  cubic  feet  per 
second  at  about  2:00  p.m.  of  June  5. 

When  analyzing  various  breach  sizes  (partial  failures),  an  elevation 
versus  discharge  curve  was  developed  for  the  downstream  boundary  condition 
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(at  the  damsite).  For  the  breach  sizes  shown  in  Plates  7  and  8  the  discharge 
corresponding  to  various  critical  depths  was  calculated  according  to1 


where 


Q  = 


(b  +  z  DC)J 
(b  +  2z  Dc) 


b  =  bottom  width  (in  feet) 

z  =  slope  of  the  sides,  horizontal  divided  by  vertical 
Dc  =  depth  (in  feet) 

g  =  acceleration  due  to  gravity  (in  feet  per  second  per  second) 
Q  =  discharge  rate  (in  cubic  feet  per  second) 


After  converting  the  various  depths  to  elevations,  the  discharge  rating 
curve  corresponding  to  each  breach  size  was  input  into  the  unsteady  flow 
program.  The  computed  discharge  hydrographs  are  shown  in  Plate  9  and 
compared  with  the  simulated  hydrograph  for  observed  conditions. 

Preliminary  information  received  indicated  that  40  percent  of  the  dam 
embankment  was  lost  such  that  an  initial  estimate  of  the  breach  configura¬ 
tion  (based  on  available  photographs)  is  shown  in  Plate  7.  The  resulting 
discharge  hydrograph,  Plate  9,  represents  an  instantaneous  failure  of  the 
size  depicted.  (The  same  is  true  in  the  case  of  the  triangular  breach  shown 
in  Plate  8.)  In  the  case  of  the  discharge  hydrograph  computed  for  observed 
conditions,  the  downstream  boundary  condition  corresponds  to  a  gradual  failure 
and  results  in  significant  differences  in  peak  and  in  the  time  to  peak. 

In  the  unsteady  flow  analysis  of  the  reservoir  system,  only  the  results 
for  n  values  of  0.04  are  included.  When  n  values  were  varied,  there  was 
no  significant  difference  in  the  computed  discharge  hydrographs  at  the 
damsite. 

Channel  Model 


In  the  case  of  the  channel,  tables  of  geometric  elements  were  prepared 
at  37  uniformly  spaced  computation  points  along  the  study  reach  (approxi¬ 
mately  4,500  feet  apart)  and  input  into  the  unsteady  flow  program.  The 


handbook  of  Hydraulics,  King  and  Brater,  Fifth  Ed.,  McGraw-Hill 
Book  Co. ,  1963,  pp.  8-11 . 
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upstream  boundary  conditions  (at  the  damsite)  are  the  various  discharge 
hydrographs  computed  in  the  unsteady  flow  analysis  of  the  reservoir. 

The  downstream  boundary  conditions  are  discharge  rating  curves  based  on 
normal  depth  calculations.  The  depth  and  corresponding  discharge  was 
determined  using  Manning's  equation 

Q  -  -±!L_  A  R2'3  S1'2 
n 


where 

n  =  coefficient  of  roughness 
A  =  area  (in  square  feet) 

R  =  hydraulic  radius  (in  feet) 

S  =  slope  (in  feet  per  foot) 

2/3 

The  values  for  the  geometric  elements,  A  and  R  ,  for  specified  depths 
(elevations),  were  taken  from  GEDA  output  for  the  most  downstream  channel 
cross  section.  The  slope  was  determined  from  available  topographic  maps. 
Values  of  n  ranged  from  0.03  to  0.07,  and  for  each  n  value  a  discharge 
rating  curve  was  developed  as  a  downstream  boundary  condition. 

In  order  to  simulate  the  actual  flow  existing  in  the  probable  flooded 
area  downstream  of  the  dam,  the  local  inflow  option  of  the  unsteady  flow 
program  was  utilized.  Major  inflow  occurred  where  Henrys  Fork  and  the 
Snake  River  entered  the  system.  The  estimated  flows  just  prior  to  failure 
were  3,500  and  5,600  cubic  feet  per  second  respectively. 

In  order  to  establish  stable  base  flow  conditions  for  each  n  value 
prescribed,  the  model  was  run  for  a  period  of  24  hours  with  the  estimated 
inflows  at  the  upstream  boundary  and  local  inflow  points  which  existed  just 
prior  to  failure.  (The  discharge  and  initial  approximations  of  water 
surface  elevation  at  each  node  or  computation  point  were  input  prior  to 
running.  Since  the  discharge  just  prior  to  failure  was  estimated  from 
observed  conditions  within  the  study  reach,  values  were  fixed,  whereas 
initial  water  surface  elevations  were  manipulated  until  stable  conditions 
resulted.)  Once  stable  base  flow  conditions  were  obtained,  the  various 
boundary  conditions  (upstream  and  downstream)  were  input  into  the  model. 

The  model  was  then  allowed  to  run  until  the  peak  of  the  flood  hydrograph 
passed  the  downstream  limit  of  study. 

Results  of  the  various  runs  are  tabulated  in  Tables  3,  4,  5  and  6. 

Plate  10  is  a  Xerox  reduction  of  the  flooded  area  resulting  from  the  June  5 
dam  failure  as  published  by  the  U.S.  Geological  Survey.  The  relative 
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TABLE  4 


COMPUTED  AND  OBSERVED  TIMES3 
FOR  LEADING  EDGE  OF  FLOOD  WAVE 


Location 

Milesb 
downstream 
from  dam 

Observed 

Arrival 

Time 

(hours) 

0.03 

Computed  arrival 
n-values 
0.04  0.05 

(hours) 

■3519 

0.07 

At  damsite 

0.0 

— 

— 

— 

— 

— 

— 

In  Teton  Canyon 

3.0 

1205 

1206 

1207 

1208 

1209 

1210 

At  Teton 

8.7 

1230 

1230 

1230 

1230 

1240 

1240 

At  Sugar  City 

12.9 

1300 

1340 

1400 

1420 

1430 

1440 

At  Rexburg 

15.8 

1340 

1440 

1500 

1520 

1530 

1540 

Henrys  Fork 
near  U.S.G.S. 
gaging  station, 
Henrys  Fork  near 
Rexburg 

21.1 

1530 

1540 

1550 

1650 

1710 

1740 

Menan  Bridge  on 
the  Snake  River 
immediately 
upstream  of  the 
Union  Pacific 

27.6 

1800 

1730 

1750 

1830 

1920 

2030 

aTimes  occurred  on  June  5,  1976. 

bDi$tances  based  on  the  unsteady  flow  data  model. 
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RESULTS  OF  UNSTEADY  FLOW  ANALYSIS 
USING  SIMULATED  DISCHARGE  HYDROGRAPH 
FOR  40  PERCENT  BREACH  SIZE 
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Elevation  table  modified  within  the  geometric  model  to  accommodate  the  higher  water  surface  elevations  resulting 
from  n  values  of  0.06  and  0.07. 
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Elevation  table  modified  within  the  geometric  model  to  accommodate  the  higher  water  surface  elevations  resulting 
from  n  values  of  0.06  and  0.07. 


location  of  the  odd  numbered  computation  points  and  conveyance  limits 
are  shown  (a  discussion  of  the  conveyance  limits  is  contained  in  the 
users  manual  for  the  GEDA  program).  (It  is  apparent  from  Plate  10 
that  the  majority  of  high  water  mark  observations  are  along  the  edges 
of  the  flooded  area.)  Maximum  water  surface  elevations  are  plotted  in 
Plates  11 ,  12  and  13. 

The  observed  high  water  marks  (taken  from  Plate  10)  were  plotted  in 
Plates  11,  12  and  13  and  designated  as  left  or  right  bank  (looking  down¬ 
stream).  The  ground  profile  shown  reflects  the  minimum  elevation  in  the 
various  cross  sections  used  in  the  geometric  model.  Horizontal  distances 
are  measured  along  the  most  probable  path  of  concentrated  flow  and  not 
along  any  particular  tributary  as  shown  in  Plate  10.  The  spacing  and 
location  of  the  odd  numbered  nodes  are  also  included  in  the  profile  plots. 
Maximum  water  surface  elevations  computed  in  the  simulation  of  observed 
data  were  plotted  for  the  odd  numbered  nodes  shown.  Water  surface  profiles 
were  then  drawn  for  each  n  value  used  in  the  analysis. 

If  a  number  of  cross  sections  are  drawn  on  Plate  10  normal  to  estimated 
flow  lines,  particularly  in  the  area  of  Rexburg  where  floodwaters  tended 
to  spread  laterally,  a  significant  variation  in  high  water  mark  elevations 
between  the  left  and  right  banks  of  a  given  cross  section  is  apparent.  In 
some  instances,  the  difference  is  as  much  as  30  feet.  This  situation  is 
also  apparent  in  Plates  11,  12  and  13.  In  general,  where  floodwaters  were 
relatively  confined  between  the  left  and  right  banks,  elevations  of  high 
water  marks  on  opposite  sides  are  essentially  the  same. 

A  comparison  of  the  observed  and  computed  maximum  water  surface  eleva¬ 
tions  (Plates  11,  12  and  13)  would  seem  to  indicate  the  following: 

(1)  Floodwaters  were  probably  concentrated  in  the  center  of 
the  flooded  area  between  the  U.S.G.S.  gaging  station  on 
the  Teton  River  and  the  downstream  limits  of  Sugar  City 
(see  Plate  10). 

(2)  Floodwaters  were  probably  concentrated  in  the  vicinity 
of  Rexburg. 

It  is  difficult  to  analyze  the  accuracy  of  the  computed  results  based  on 
observed  high  water  marks.  Any  detailed  analysis  would  definitely  require 
high  water  marks  within  the  limits  of  flooding.  If  such  could  be  obtained 
a  greater  effort  could  be  made  in  stipulating  conveyance  limits  and  in  the 
actual  calibration  of  an  unsteady  flow  model. 

Tables  1  and  4  compare  the  estimated  (based  on  field  observations)  and 
computed  arrival  times  of  the  leading  edge  of  the  flood  wave,  whereas 
Tables  2  and  5  compare  estimated  (based  on  field  observations)  and  computed 
maximum  discharges  at  particular  locations  within  the  study  area.  Plate  14 
gives  an  indication  of  the  time  required  for  passage  of  the  flood  wave 
(based  on  n=0.04). 
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In  the  unsteady  flow  model,  determination  of  the  arrival  time  of  the 
leading  edge  of  the  flood  wave  depends  on  a  noticeable  difference  in  the 
computed  water  surface  elevation  occurring  at  a  given  node.  Since  minor 
fluctuations  of  the  water  surface  were  present  in  the  base  flow  and  since 
the  printout  interval  was  10  minutes  after  12:30  p.m. ,  arrival  times  were 
estimated  (probably  ±  30  minutes). 

Computed  peak  discharges  (Table  5)  differ  significantly  from  those 
estimated  from  observed  data  (Table  2).  It  is  reasonable  to  assume  though 
that  the  discharges  of  Table  2  are  only  rough  approximations. 

In  the  development  of  the  unsteady  flow  data  model,  there  was  no  effort 
made  to  locate  computation  points  at  those  locations  specified  in  Tables  1 
and  2.  (Usually,  computation  points  are  located  such  that  results  are  printed 
at  desired  locations.) 

5.  STEADY  FLOW  MODEL 

The  information  originally  coded  in  the  HEC-2  format  and  utilized  in 
development  of  the  channel  geometric  model,  was  used  to  determine  the  storage- 
outflow  relationships  for  various  reaches  downstream  of  the  dam.  Multiple 
profiles  were  run  using  HEC-2  and  the  volumes  in  a  given  reach,  corresponding 
to  a  specified  discharge,  computed.  (It  should  be  noted  that  an  n-value  of 
0.04  was  used  in  the  multiple  profile  computations.) 

The  discharge  hydrograph  (see  Plate  6)  based  on  observed  data,  was  routed 
downstream  with  the  storage-outflow  relationships  previously  determined. 

Channel  routing  was  accomplished  with  HEC-11  using  modified  Puls.  The  peak 
discharges  computed  for  each  reach  were  subsequently  input  into  the  HEC-2 
data  deck. 

Before  running,  the  HEC-2  deck  was  modified  to  reflect  the  conveyance 
limits  adopted  in  the  unsteady  flow  analysis  (see  Plate  10).  Encroachment 
limits  were  set  with  the  X3  card  as  close  as  possible  to  the  conveyance 
limits  previously  designated  (this  was  necessary  since  GR  stations  did  not 
always  correspond  to  the  conveyance  limits  specified  in  the  geometric  data 
model).  The  resulting  water  surface  profile  is  shown  in  Plates  15,  16  and  17 
(superimposed  for  comparison). 

The  difference  in  the  computed  water  surface  profiles  between  the  two 
methods  (the  steady  flow  and  the  unsteady  flow  analysis)  is  due,  in  part,  to 
the  difference  in  the  computed  discharges.  Discharges  computed  using  modified 
Puls  were  lower  in  the  upstream  channel  reaches  and  higher  in  the  lower 
reaches.  This  is  attributed  to  the  significant  effect  of  the  inertia  component 
(in  the  unsteady  model)  inmediately  downstream  of  the  dam  and  in  the  way 
storage  is  treated  in  the  two  methods. 

1 HEC-1  Flood  Hydroqraph  Package,  Generalized  Computer  Program  (723-010), 
Hydrologic  Engineering  Center,  U.S.  Army  Corps  of  Engineers,  January  1973. 
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APPENDIX  B 


Guidelines  for  Analyzing  a  Dam  Break  Flood  with  the  Computer  Program 
"Gradually  Varied  Unsteady  Flow  Profiles" 


1 .  Data  Requirements 

a.  Topographic  maps  are  required.  The  area  of  coveraqe  should  extend 
from  the  upper  end  of  the  reservoir  in  question  to  the  most  downstream  point 
of  interest.  The  contour  interval  should  be  less  than  the  level  of  accuracy 
expected  from  the  study. 

b.  Aerial  or  ground  photographs  are  desirable.  These  provide  hydraulic 
engineers  with  a  basis  for  estimating  hydraulic  roughness  coefficients. 
Consequently,  they  are  not  needed  for  the  reservoir  area  but  for  the  channel 
and  flood  plain. 

c.  Reservoir  capacity  curves  or  tables  are  required  and  may  be 
calculated  from  the  topographic  maps.  These  show  the  volume  of  water  stored 
as  reservoir  depth  increases. 

d.  Breach  size  and  shape  is  required  and  may  be  estimated. 

e.  Reservoir  level  at  breaching  is  required  and  may  be  estimated  from 
aerial  photography  if  reservoir  records  are  not  available. 

2.  Assumptions. 

a.  It  is  customary  to  assume  that  the  breach  is  developed  instantly 

and  completely  to  the  desired  size  and  shape.  Otherwise,  several  Instantaneous 
breaches  of  different  sizes  may  be  postulated  and  the  resulting  floods 
calculated. 

b.  Critical  depth  controls  at  all  partial  breaches.  In  view  of  the 
values  presented  in  Table  2,  main  body  of  this  report,  inertia  is  not  a 
major  consideration  over  a  wide  range  of  partial  breach  sizes. 

c.  The  streambed  and  banks  are  fixed  (i.e.,  do  not  erode)  during 
the  event. 

d.  Any  bridge  across  the  stream  fails  instantly  upon  impact  of  the 
flood  wave.  The  resulting  energy  loss  is  negligible. 


3.  General  Procedure.  The  general  procedure  depends  somewhat  upon  the 
type  of  breach  postulated.  That  is,  if  the  entire  structure  is  removed, 
the  final  step  of  the  analysis  will  calculate  the  outflow  hydroqraph  from 
the  reservoir  plus  route  that  outflow  hydroqraph  to  all  downstream  points 
of  interest  in  a  sinqle  operation. 

If  only  a  partial  breach  is  postulated,  the  outflow  hydrooraph  is  calculated 
first  and  the  routinq  to  downstream  points  is  accomplished  in  a  separate 
operation.  The  same  analytic  technique  is  used  for  both. 

The  obvious  advantaqe  of  the  first  approach  is  a  savinqs  in  the  time 
required  for  performing  the  analysis.  A  mor*»  siqnificant  advantaqe  is 
that  the  method  includes  any  submeroence  effects  due  to  the  rise  in 
tailwater  elevation  at  the  dam.  This  submergence  effect  can  be  suffi¬ 
ciently  large  to  change  the  shape  of  the  outflow  hydroqraph  from  the 
reservoir. 

partial  breaches  cannot  be  analyzed  using  the  first  approach  because  the 
existinq  computer  program  will  not  accommodate  a  hydraulic  jump.  On  the 
other  hand,  the  rise  in  tailwater  elevation  is  usually  small  enough  so 
that  no  submerqence  effect  develops.  Therefore,  the  second  atmroach  is 
adequate. 

4.  Calculation  of  Outflow  Hydroqraph  Using  the  Partial  Breach  Approach. 

The  details  which  follow  present  the  partial  breach  approach.  The  first 
step  in  the  calculation  will  determine  the  outflow  hydroqraph,  paragraphs 
a-f.  Routing  the  flood  hydrograph  to  downstream  points  is  presented  in 
paragraph  8. 

a.  Calibration  of  Cross  Sectional  Data  to  Reservoir  Volume.  Locate 
cross  sections  on  the  topographic  map  so  reservoi r”volume  may  be  calculated 
using  the  average  end  area  technique.  Include  all  major  tributary  arms  in 
the  reservoir.  Code  this  data  for  the  computer  program  "Geometric  Elements 
from  Cross  Section  Co-ordinates"  (GEDA)  and  calculate  the  reservoir  volume 
for  the  range  of  elevations  up  to  full  pool.  This  result  should  agree  with 
the  reservoir  capacity  curve.  If  It  differs,  (  +  5"),  cross  sections  which 
are  more  representative  should  be  developed.  Usually,  this  means  recoding 
some  of  the  tributary  arms. 

b.  Geometric  Model  for  Reservoir.  Select  a  distance  between  computation 
nodes  that  comnlies  with  criteria  in  the  routing  program,  "Gradually  Varied 
Unsteady  Flow  Profiles"  (USTrLO).  Values  typically  range  between  3- mile 

(408  m)  and  5  miles  (15,000  m).  Execute  GEDA  to  produce  the  required  tables 
of  geometric  properties  and  n-values  at  each  computation  node.  The  most 
downstream  section  should  be  at.  the  dam  axis  and  include  the  breach  that 
was  postulated. 
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c.  Downstream  Boundary.  Calculate  a  critical  depth  ration  curv^ 
for  the  breach  size  and  shape  postulated.  Each  elevation  selected  will 
establish  an  area  and  top  width  for  use  in  the  followinq  equation. 

0  =  A  -ylq  0 
cr  cr  ’  cr 


where: 

A  „  =  cross  sectional  area,  square  feet 
cr 


D 


cr 


A  /R 
cr 


cr 


R  „  =  too  width,  feet 
cr 


o  =  acceleration  of  qravitv  feet/socond/second 


d.  Initial  Conditions.  Establish  the  initial  water  surface  elevation 
and  discharqe  at  each  computation  node  in  the  reservoir.  It  is  exoedient  to 
always  select  a  horizontal  water  surface  and  zero  discharqe.  Although,  one 
may  oostulate  whatever  conditions  he  desires  bv  startinq  the  comouter  model 
at  the  above  conditions  and  simulatinq  inflow/ outflow  records  ud  to  the  time 
when  breaching  is  anticipated. 


e.  Upstream  Poundary.  Usually,  the  initial  reservoir  inflow  is 
assumed  to  be  zero.  It  may  be  otherwise  if  so  desired.  In  either  case, 
code  thp  inflow  hydroqraph  as  described  in  the  USTFLO  users  manual. 


f.  Calculating  the  Outflow  Hydroqraph.  Select  an  interval  between 
printouts  that  is  very  short  (a  minute  or  less)  during  the  early  part  of 
the  hydroqraph  (5  minutes  or  so).  Usually,  the  calculated  peak  will  occur 
during  this  interval  and  the  time  between  printouts  may  be  increased  to  1° 
or  15  minutes.  In  any  case,  obtain  a  sufficient  amount  of  Printout  to 
define  the  entire  outflow  hydroqranh  shape. 


5.  Routin')  th_e_  bam  RreaF  Flood 

a.  Geometric  'icdel  fo_r  Houtjni  Downstream.  Locate  cross  sections  to 
lhp  downstream  points  of  interest  us  inn  information  oained  from  the  reservoir 
model  and  followinq  s‘eps  in  the  users  manual.  Select  the  distance  between 
computation  nodes.  (This  distance  does  not  have  to  be  the  same  as  that  usM 
in  the  reservoir  when  routin')  with  a  two  step  process.  However,  use  the 
sane  qeneral  guidelines.)  Code  tbe  data  and  execute  CFHA  to  produce  the  re¬ 
quired  geometric  data  tables. 
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b.  Downstream  Boundary.  If  possible,  set  the  downstream  boundary 
a  couDle  of  nodes  beyond  the  ooint  of  interest.  Calculate  a  rating  curve 
by  slope-area  or  by  extrapolating  from  observed  values. 

c.  Upstream  Boundary.  The  outflow  hydrograph  calculated  in  paragraph 
4  becomes  the  upstream  boundary  for  this  step. 

d.  Initial  Conditions.  Establish  initial  conditions  by  calculating 
a  base  flow  profile.  (?ero  or  negative  depth  is  not  an  acceptable  initial 
condition  .  )  It  is  expedient  to  prescribe  a  bank  full  water  discharge,  to 
let  the  computer  model  stabilize  by  simulating  the  steady  flow  profile  for 
that  discharge  and  finally  to  decrease  the  inflow,  gradually,  until  the 
desired  base  flow  discharge  is  reached.  Simulate  the  base  flow  conditions 
for  a  sufficiently  lonq  period  of  time  to  establish  a  steady  flow  profile. 
Code  this  profile  (water  surface  elevation,  discharge)  to  form  the  initial 
conditions  for  routing  the  dam  break  flood. 

e.  noutinfi  the  Flood.  The  computer  orogram  referenced  in  subparagraph 
lc  produces  both  a  discharge  hydrograph  and  a  water  surface  elevation  hydro¬ 
graph  at  every  computation  nod^.  A  summary  table  gives  maximum  and  minimum 
values  for  elevation,  discharge  and  velocity. 
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